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We consider a version of the McLerran-Venugopalan model by Lam and Mahlon 



O' 

CSl I where confinement is implemented via colored noise in the infrared. This model does 
> ■ 

O ■ not assume an infinite momentum frame nor that the boosted nuclei are infinitely 

O: 

. thin; rather, nuclei have a finite extension in the longitudinal direction and therefore 

" depend on the longitudinal coordinate. In this fully three dimensional framework 

■ an X dependence of the gluon distribution function emerges naturally. In order to 

■ fix the parameters of the model, we calculate the gluon distribution function and 
r~| . compare it with the JR09 parametrization of the data. We explore the parameter 



space of the model to attain a working framework that can be used to calculate the 
initial conditions in heavy ion collisions. 
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INTRODUCTION 

In heavy ion collisions at very high energies, soft (small-x) partons, mostly gluons, consti- 
tute a big portion of the wave function of the colliding nuclei. Since the occupation number 
of the gluons is very high, they can be treated in the framework of classical Yang-Mills 
theory. When two nuclei pass through each other, the classical nonabelian fields from each 
nucleus interact with each other and form color flux tubes between the receding nuclei, which 
ultimately lead to the quark gluon plasma. The energy density distribution of the initial 
state determines the multiplicity and spectrum of the final particles seen in the detectors. 
The initial energy density distribution is an input to the hydrodynamics description of the 
evolving quark gluon plasma and, in principle, it can be calculated ab initio from Quantum 
Chromodynamics (QCD). 

The energy initially deposited in the color flux tubes can be related to the two-point vector 
potential correlation function {Af{x)A'i^{x')) where the setup may be two (transverse) or 
three dimensional, depending on whether or not the Lorentz contracted nucleus is assumed 
to have a longitudinal thickness in the lab frame. This correlator can be derived analytically 
from the color charge density correlator {p"-{x)p^{x')) by using the classical Yang-Mills 
equations. 

The color charge density of a nucleus p°'{x) during a collision cannot be known; on 
the other hand, the fluctuations in the color charge density can be studied in the effective 
field theory approach with ensemble averaging. In the Color Glass Condensate (CGC), a 
framework for slowly evolving high density gluons within the ultrarelativistic nucleus, the 
fast partons are seen as sources of the small-x (soft) gluon radiation, where x is the fraction of 
the total longitudinal momentum carried by a parton. After integrating out the fast partons, 
the observables can be calculated by averaging them over the ensemble of all possible color 
charge configurations. Once the correlator {p°'[x)p^[x')) is specified, it can be linked to the 
vector field correlator {A'^{x)A\{x')) from which the gluon distribution function xg{x,Q'^) 
and other observables, such as the initial energy density, can be calculated. Here (■■■) 



denotes the ensemble average. See Ref. 



11] for a comprehensive review. 



In the early formulation of the CGC by McLerran and Venugopalan , the spectrum 
of the Gaussian color fluctuations was taken to be white noise. This results in arbitrarily 
long wavelength fluctuations where {A'^[x)A'l[x')) diverges in the infrared. This problem 
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originates from the fact that confinement effects (color neutrahty) have not been taken into 
account. Adding a gluon mass to the gluon propagator so as to bypass this problem breaks 
gauge invariance which is required later to convert the solution of the classical Yang-Mills 
equations from the axial gauge to the light-cone gauge by means of Wilson lines. 

In the original McLerran-Venugopalan (MV) model, the colliding nuclei are considered 
to be two-dimensional infinitely thin sheets traveling at the speed of light. In this ap- 
proximation, the model does not depend on the longitidunal coordinate, hence there is no 
dependence on the momentum fraction x of a given parton. In other words, all of the par- 
tons in the nucleus have the same momentum fraction x. Being an artifact of the infinite 
momentum frame, the lack of x dependence in the model does not reflect the true nature of 
the X dependent gluon distribution functions. 

In this paper, we consider a color neutral three dimensional McLerran-Venugopalan 

n h 

(SdMVn) model by Lam and Mahlon [7|, 18|, where the spectrum of the Gaussian fluctu- 
ations is taken to be infrared-safe colored noise. Colored noise creates a correlation hole in 
the correlation function on the size of a nucleon and this leads to color neutrality. (The term 
"colored noise" is not related to the color charge of QCD.) This model comes with another 
important feature that the incoming nuclei are not exactly on the light cone. Longitudinal 
coordinate dependence produces an x dependent gluon distribution function. This enables 
us to compare the model with the data parametrization of gluon distribution functions and 
therefore fix the parameters of the SdMVn model. Once this is achieved, the SdMVn model 
can be used to calculate the initial energy density distribution in heavy ion collisions. 

The MV model provides a framework for calculating the vector field correlation function 
{A'^{x)A'^{x')) from which the gluon distribution function as well as the initial energy density 
distribution can be calculated. We emphasize that our ultimate goal is calculating the latter, 
which will be pursued in a follow-up paper. In this paper we focus on the calculation of 
gluon distribution function and its comparison with data for the purpose of fixing the free 
parameters in the SdMVn model. 

The outline of the paper is as follows. First, we give a brief overview of the color neutral x 
dependent McLerran-Venugopalan (SdMVn) model given by Lam and Mahlon [3,0]. Next, 
we calculate the gluon distribution function from the SdMVn model. This model works at 
low-Q^ where no data for the gluon distribution functions is available. For that reason, we 
evolve our results with the DGLAP equation to higher momenta where data is available. 
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This is followed by a discussion of the data parametrization that we use to compare with the 
3dMVn model. Lastly, we compare the 3dMVn with the data parametrization. The final 
section summarizes our analysis and includes a discussion regarding how the results will be 
used in the future to calculate initial energy density distribution. 



COLOR NEUTRALITY 



In the original two dimensional MV model, the average and the fluctuations of the color 
charge density of a nucleus are given by 

(P"(^x)) = 0, (1) 
{p'^{x^)p\x'^)) = S^'plVix^- x'^), (2) 

where Xj_ is the transverse coordinate system for the infinitely thin nucleus, p\ is the average 
color charge density squared per unit (transverse) area for a nucleus, and (■ ■ ■ ) denotes the 
ensemble average. The function V determines the spectrum of the fluctuations. When white 
noise is assumed, T)( ) = — x'j_), no correlation occurs between different points 

in a nucleus. This also means that fluctuations in each momentum mode, including the zero 
mode, are equally likely and there is no characteristic scale. In this case, even though the 
correlator {A^{xj_)A'l{x'j_)) is calculated with an infrared cutoff Aqcd Iq], it still diverges like 
(x^)^j- in the infrared. Therefore, the Fourier transform of it, which is necessary to calculate 
the gluon distribution function, does not exist p'J. This cut-off can also be understood as 
a gluon mass in the gluon propagator in the form {q^ + "^gluon)~^• explicitly breaks the 
gauge invariance that is needed to convert the solutions of classical Yang-Mills equations 
from the axial gauge to the light-cone gauge where Wilson lines are to be used. 

In reality, a nucleus is color neutral on scales larger than a nucleon size. A color neutral 
correlation function can be derived from a simple model of a nucleus where nucleons are 
composed of quark and antiquark pairs [l^. If we consider a two dimensional nucleus for 
a moment, the assumption that only the quark and antiquark pair from the same nucleon 
can interact with each other produces a correlator of the form 7| 

{p^{x^)p\x'^)) = 5^'p\ Y^""^ - <^ - • (3) 

Here the parameter A, which we will refer to as a correlation length, is of the order of Aq^j-,. 
This parameter takes away the need for a sharp infrared cutoff as it is used in the original 
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MV treatment; moreover, it makes the Fourier transform of the vector field correlation 



function well-defined and it renders the model infrared safe. In Eq. ([3]), decorrelation due 
to the white noise T)( = 6'^{xj_ — x'j_) is modified by the last term for distances 

I < A ~ ^QCD' reflecting the assumption that there is a correlation between the 
partons confined to a region smaller than the nucleon size. The second part in Eq. ([3]) 
removes the zero mode \q^\ = from the white noise spectrum, hence colored noise. The 
Fourier transform of the correlator in Eq. (|3]) is proportional to 1 — e"^^*^?- , which vanishes 
as |q_|_| — > 0. This makes the model infrared safe. 

In this paper, we adopt a three dimensional correlator given by Lam and Mahlon js] 

3 exp{-^/3\x\/X) 



(p'^(0, 0)^(0^)) =5»''4 



47rA2 



\x\ 



(4) 



where k\ is the three dimensional average color charge density 



3ACp 1 



3A 1 



(5) 



Here is the number of colors, V is the volume of a nucleus and the color factor is defined 
as Cp = {N^ — l)/{2Nc). We define d^x = dx\\d'^Xj_ where x\\ is the longitudinal coordinate 
in the direction of the beam axis, and Xj_ is the coordinate on the transverse plane. The 
Fourier transform of the corralator (j4]) is given by 

1 



1 



(6) 



1 + A2gV3_ 

Despite the slight difference between the correlators in Eqs. (E]) and (jl]), they produce 
similar results. It can be readily seen from Eq. 
spectrum since (p"p^) — t- as \q\ — )■ 0. 



that the zero mode does not exist in the 



LONGITUDINAL DEPENDENCE 



In ultrarelativistic heavy ion collisions, the infinite momentum frame provides a good 
starting point for the parton picture of the nucleus. However, the nucleus becomes infinitely 
thin when it is exactly on the light cone and the gluon distribution function turns out to 
be independent of x. Here we review the formulation by Lam and Mahlon [8] for the case 
where the nucleus is boosted to speed /3. In the lab frame, the thickness a nucleus of radius 
R becomes of the order of R/j. 
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The current for a nucleus moving in the +z direction is 

J° = p(-^,a;^J; J, = 0, (7) 

where the subscript "r" stands for the rest frame and the negative sign in front of the z is 
for later convenience. With the redefinition = (t ± z) / -\/2, the current can be 

rewritten in the light-cone coordinates (still in the rest frame), 

When we go from the rest frame to the lab frame where the nucleus moves with speed (3, 
the current in Eq. ([8]) becomes 

J+ = -p ( -X" - xA; .r = J J+; = 0, (9) 



where 



Here we can define a new longitudinal coordinate 

X\\ = -x~ — -x~^, (11) 

which is essentially the Lorentz tranformation of —z = {x~ — x~^) / \/2 (see the definition 
below Eq. ((Tj)). The Eq. can be contrasted with the commonly used current in the 
infinite momentum frame (imf) where the nucleus is infinitely thin, 

J+f = 5(x-)p(a;J; Jr^, = 0; Jximf = 0. (12) 



CALCULATION AND EVOLUTION OF THE GLUON DISTRIBUTION 

FUNCTION 

We now turn to the calculation of the gluon distribution function from the 3dMVn model. 
Later we will compare it with data in order to fix the parameters and A. The gluon 
distribution function of a nucleus with a baryon number A can be expressed as an integral 
of the gluon number density over the transverse momenta. 
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(14) 



The gluon number density in the lab frame for a nucleus moving with speed /3 is given as a 

Fourier transform of the two-point vector field correlation function 

dN _ q\\ 
dq\\d'^q± 4:n^ 

Here q\\ is the momentum conjugate to the coordinate defined in Eq. ( ITTl) . The Eq. iHM 
can be related to dN / dxd^q^ by using the relation x = qw/m, which itself can be derived 
from the definition x = q'^ /Q^ = eq'^/m. Here m is the nucleon mass, g"*" and Q'^ are the 
momenta of the gluon and the nucleon. 

Solving the classical Yang-Mills equations for the source in Eq. ([8]) in the covariant 
gauge and transforming the solution into the light-cone gauge by using Wilson lines, one 
can express the correlator {A^{x)A'^{x')) in Eq. (HM in terms of the color charge density 
correlator {p°'{x)p^{x')). We skip the details of this calculation here and refer the reader 
to the original paper Using the correlator in Eq. (jl]), the gluon number density can be 
written as 

= / rf^A^e--^-£(x; A^)S{v'LiA^)), (15) 

Here the integration is over the transverse coordinate A_l The functions C 

and L are convolutions of the two gluon propagators at two different points in the same 
nucleus connected by the three dimensional noise term V{x — x'). These two functions can 
be seen as a pair distribution function. The nuclear correction factor S takes into account 
the nuclear geometry and f ^ controls the strength of dependence on this geometry. 



The functions used in Eq. (1151) are given as 

1 



C{x; Aj 



and 



L(A,) = — 



A_2 

677 



Kr 



127r' 
f 73 A 



8( 

xmA)^-K'o(xmA_L), 



+ In 



( V3A_ 
2A 



+ 7e 



(16) 



(17) 



In the calculations of Lam and Mahlon [85], the nuclear matter density is taken to be uniform. 
The nuclear correction factor S is determined by geometry of the nucleus. For cylindrical 
and spherical nuclei it is given by 

1, 



S(z) 



1) 



(cylindrical) , 



(18) 



-[1 - -z^ + e''(z - 1)] (spherical). 
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This function appears in Eq. (fT5|) as £^(f ^L(A_|_)). The controls how much the results 
depend on the nuclear geometry; it is given by 

^ ^ 247ra2Ai/3A^(,j3 (cylindrical), 



(19) 

^ ^ 367ca^,A^/^AlcD (spherical) . 

The difference between the spherical and cylindrical nuclei is negligibly small. In our calcu- 
lations, we use the formulae for a cylindrical nucleus. 

Now we are ready to calculate the gluon distribution function xg^lXjQ"^) for a nucleus 
by substituting Eq. ( ITSl) in Eq. ( fT3l) . 



xgAix, Q^) = ^ / <fq^ [ rf2A^e"^^-^^£(x; A l.)£{v^L{A^)). (20) 

J\q^\i:Q J 

The momentum integration in the equation above is trivial. Owing to the longitudinal 
coordinate x\\, the fractional momentum distribution function xgAiXjQ"^) in Eq. fl20|) is x 
dependent, in contrast to the results of the original two dimensional MV model. 

The region of validity of the 3dMVn model within the two dimensional parameter space 
(x, Q^), likewise the original MV model, is restricted by the assumptions made regarding 
the weak coupling limit, color coherence and color averaging js]. 

The lower limit of x is set for a gold nucleus by 

^-1/3 

X ~ 0.035 (21) 

ma 

where a is the nucleon size ~ 1 fm. This is the longitudinal momentum fraction at which 
the gluons start to resolve the Lorentz contracted thickness of the nucleus. 
The weak coupling limit imposes an upper limit on x such that 

4 



< 



0.25. (22) 



TTma 

For very large Q^, not enough color charge would be seen, hence the Gaussian approximation 
for the color charge fluctuations {p"'{x)p^{x')) would not be valid. Similarly, at the scale of 
color neutral nucleons where is very small, again there would not be enough color charge 
to average over. These considerations put limits on the momentum scale 

j^x <Q'< -^i. (23) 

4 ma X 
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X 



FIG. 1. Approximate region of validity of the 3dMVn model for gold nucleus (A=197). 
For Aqcd = 200 MeV, it becomes 

0.16x<g^<^. (24) 

X 

In our calculations, we will take = 0.55 GeV^. The reason for this choice is that it 
is at the upper limit of validity of the 3dMVn model and at the lower limit of the parton 
distribution function to which we will compare. While the model spans a larger range in x 
for Q'^ < 0.55 GeV^ (see Figured]), no data is available at those scales. At the energy scale 

= 0.55 GeV^, x is restricted to the range 

0.035 <x< 0.060. (25) 

Figure [1] shows the approximate validity region of the SdMVn model for a gold nucleus. We 
consider gold because it provides for the greatest range of validity of the model as represented 
by Eq. fl^T]) . The results for lead would be indistinguishable. 

The parameters of the SdMVn model need to be fixed before performing the numerical 
integration in Eq. (120]) . We take the nucleon mass m = 0.94 GeV, Aqcd = 0.2 GeV and A = 
197. The coupling constant as and the correlation length A are treated as free parameters. 
We evaluate xgA{x,Q^) numerically at = 0.55 GeV^ in the range 0.035 < x < 0.060 for 
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several values of and A in the range 



0.1 <as < 1, 

- - (26) 
0.2 fm <A < 2.6 fm. 



Then we fit xqaIx, Q"^) to a form 

a;r?A(x,Q') = 6x^(1 -x)^ (27) 

and find {6, c, c?} for each set of {as,A}. After the numerical integration, we obtain 
x5'i97(x,Q^) for the whole nucleus for several points in the parameter space {a5,A}. The 
scale = 0.55 GeV^ at which we calculated xqa^x, Q^) is the lower limit of the JR09 data 
parametrization for the parton distribution functions. For this reason, we will evolve the 
gluon distribution functions evaluated at = 0.55 GeV^ to higher momentum values such 
as = 25, 100 and 1000 GeV^ and compare them with data at those scales. 

Gluon distribution functions are evolved to different energy scales with the DGLAP 



equation 
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14| . For this purpose, we employ the code QCDNUM17 [15||. We will assume 



that nucleus is a dilute system of nucleons and take the nuclear modification factor = 1 
Q. Hence, 

xgp/A{x, Q^) = ^xgA{x,Q'^), (28) 

where the subscript p/A refers to a nucleon in a nuclear environment. Henceforth, we will 
use the gluon distribution function for a single nucleon in a gold nucleus. It should be 
kept in mind that 5'p/^(x, Q"^) is in principle different than a gluon distribution function for 
an isolated proton gp{x, Q"^) since the nuclear effects introduce enhancement or shadowing 
depending on the energy scale . 

We run QCDNUM in the variable-flavor number scheme (vfns) and at the NNLO(MS) 
order. For the DGLAP evolution of the gluons, valence and sea quark distributions at 
the initial scale should also be provided to the evolution code. At the initial scale = 
0.55 GeV^, the calculated gluon distribution function is the main input to QCDNUM. As 



for the valence and sea quark distributions, the JR09 17| data parametrization is utilized. 



The evolution is repeated for several values of and A in the range given in Eq. (|26|) 
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DATA PARAMETRIZATION 



We will compare our results with the parametrization of data by JR09 17|. We use the 
pdfs that work in the vfns scheme and at the NNLO(MS) order. The JR09 is valid in the 
ranges 

0-55 < iqS GeV^, 

(29) 

10"^ <x < 1. 



The JR09 provides pdfs for separate nucleons xgp{x, Q"^) but it does not contain any 
information about the nuclear modification function R^. Our results, on the other hand, 
are for a nucleon in a nuclear environment XQp/Ai^, Q"^)] therefore, Ra should be, in principle, 



taken into account. However, the nuclear pdfs (for example nCTEQ 18[) are available only 
for > 1 GeV^. As pointed out earlier, since the DGLAP equation mixes quarks and 
gluons during the evolution, we need the valence and sea quark distribution functions at 
the initial energy scale so that the gluon distribution function calculated from the model 
can be evolved to the higher energies. Hence, for consistency, we will utilize the quark pdfs 
provided by JR09 at = 0.55 GeV^ as an input to the DGLAP evolution. The discrepancy 



between JR09 and nCTEQ 18|, encoded in the nuclear modification factor, is at most 5% 
at Q'^ = 25 GeV^. At Q'^ = 100 and 1000 GeV^ the discrepancy is negligibly small. For that 
reason, we think the lack of information regarding Ra will not cause a significant error in 
our analysis. 



BEST FIT PARAMETERS 



Now we seek the sets of parameters {a^. A} that produce the best fits. For this reason, we 
compare the gluon distribution functions, calculated and evolved to = 100 GeV^, with 
JR09 at the same energy. 

In the DGLAP evolution to higher energies, the contribution to the radiation at a par- 
ticular value of X always comes from the sources from the larger x region. Hence, the effect 
of gluons calculated from 3dMVn in the range 0.035 < a; < 0.060 at = 0.55 GeV^ will be 
more prominent at smaller x values at = 100 GeV^. For this purpose, we will compare 
the evolved model with JR09 in the range 



0.015 <x< 0.04. 



(30) 
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X 

FIG. 2. Comparison of the gluon distribution function of a nucleon from the 3dMVn model for 
various values of {as,X} with the JR09 data parametrization at = 100 GeV'^. The horizontal 
axis is in logarithmic scale. For distances larger than the correlation length A, gluons are not 
correlated and hence the color neutrality condition is satisfied. The model is reliable in the fit 
region bounded by the two vertical lines. The uncertainty in JR09 is shown with an error band 
and it is about 5%. The discrepancy between the best fit curve (a^ = 0.5 and A = 1.8 fm) and the 
JR09 in the fit region is only 2%. 

The gluons from x < 0.015 do not enter the DGLAP evolution and hence they do not 
contribute to the region in Eq. ( 130|) . The gluons from 0.015 < x < 1 will enter the DGLAP 
evolution. We argue that the parametrization in Eq. fl27p is a good approximation for the 
large-x gluons. 

Our analysis shows that the best fit occurs at the values = 0.5 and A = 1.8 fm. In 
Figure [2] we show a comparison between the JR09 data parametrization and 3dMVn model 
at = 100 GeV^ for as = 0.5 and A = 1.8 fm. We find almost identical plots at other 
values in the range 0.015 < x < 0.04. 

At other values of and A, we find that the model underestimates the data. As ag 
increases, the 3dMVn curve increases and the disagreement between the model and JR09 



13 



decreases. For values as > 0.5 the model does not change significantly. Hence, we take 
as = 0.5 since this is the smallest value of for which a good fit can be obtained. Inter- 
estingly, this "freezing" behavior of the str ong coupling constant particularly at as = 0.5, 



imposed in other approaches by hand (see 19|, |20|), arises naturally in the 3dMVn model. 



Adjusting A changes the model curve only slightly. For as = 0.4 and values of A greater 
than 1.8 fm, we do not find a fit as good as the one for = 0.5 and A = 1.8 fm. 

There are various sources of uncertainty in the MV model which are shared with the 
3dMVn model. In the paradigm of strong classical color fields, one wishes to employ solutions 
of the classical Yang-Mills equation which is given in covariant gauge by 

(Vi + d^)A'' = gr + 2^^[A^ d.A'] - ig[A,, d'' A^] + /[A^, [A^ A^]], (31) 

where 9y = — 29+(9_. The second and third terms on the right hand side of the Eq. f l3T|) 
are responsible for the three gluon interaction vertex, and the last term is for the four gluon 
interaction vertex in the quantized theory. In the framework of Color Glass Condensate, 
only the linearized version of Eq. (13T]) is used. The nonabelian feature of the model sets in 
when the solutions of the linearized equation in the covariant gauge are transformed to the 
ight-cone gauge by means of the full gauge transformation including the nonabelian term 

3 

Ai^{x) = U(x)A^f,,U-i(x) - ^[9^U(x)]U-i(x). (32) 

The reason why the nonlinear terms in Eq. flHT]) are omitted is simply because the solutions 
of the fully nonlinear theory are not known. In addition, the Green's function for gluons, 
which is essential in the calculation of {Af{x)A'l{x')) from {p"'{x)p''{x')) , can only be defined 
in the linear theory. The gluon propagator is calculated from 

iVl + dl)Gix) = 6%x). (33) 

The origin of the infrared divergences in the MV model lies here and this is also why 
confinement effects need to be introduced by hand by using colored noise. In principle, 
the correlation function with colored noise in Eq. (jl]) should be calculable from the fully 
nonlinear theory. 

The linear approximation can only be justified if ^ 1 so that the nonlinear terms can 
be neglected. However, the strong couphng constant is not expected to be so small and 
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nonlinear terms are as important as the source term. Hence, the hnear approximation of 
Eq. fl3Tl) is the main source of the uncertainty in any reahstic analysis. 

Another source of uncertainty may be the assumption that the nucleus is taken to be 
much larger than a nucleon, R ^ Aq^j-,. For a gold nucleus, the corrections may be as large 
as 17% since Aq^j-,//? ~ A~^^^ ~ 0.17. Also quantum corrections at smaller values of x may 
be important even in the weak coupling limit. These corrections are of the form a,ln(l/a;) 

n n n 

and discussed in the references [5|, 16|, 18[. Lastly, although we estimate that the error due to 
neglecting the nuclear effects should be small, it would be worth of repeating the analysis 
presented in this paper by using nuclear pdf 's once they are available for ~ 0.55 GeV^. 



SUMMARY AND OUTLOOK 

The framework of Color Glass Condensate allows one to calculate the vector field cor- 
relation function. It can be used to calculate the gluon distribution function of the ultra- 
relativistic nuclei and the initial energy density distribution due to the interacting classical 
color fields produced by the colliding nuclei in heavy ion collisions. 

We have examined a three dimensional color neutral version of the McLerran-Venugopalan 
model. The 3dMVn model is finite in the infrared and therefore an infrared cutoff is not 
needed. In addition, the results of this model are x dependent due to the intrinsic three 
dimensional treatment of the nucleus in contrast to the approximation of infinitely thin 
nucleus. 

In order to explore the parameter space of these two variables, we have calculated the 
gluon distribution function for several values of as and A. Then we compared our calculations 
with the JR09 parametrization. We have found the best fit between the gluon distribution 
function from the JR09 parametrization and the one calculated from the 3dMVn model 
occurs at = 0.5 and A = 1.8 fm. At other values of these two parameters the model 
underestimates the data. This may be due to the uncertainty in the assumptions and rough 
estimates made during the construction of the model as well as the other uncertainties 
discussed in the previous section. 

In this work, the ultimate goal was not calculating the gluon distribution function per se, 
but to attain a working model for further use, particularly in the calculation of the initial 
energy density distribution in heavy ion collisions. This work is in progress. In the future 
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work, we plan to calculate the initial energy density in heavy ion collision in the framework 
of the 3dMVn model. 
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